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There is a reformulation of magnetohydrodynamics in which the fundamental dynamical quan- 
tities are the positions and velocities of the lines of magnetic flux in the plasma, which turn out to 
obey equations of motion very much like ideal strings. We use this approach to study the evolution 

- - - of a primordial magnetic field generated during the radiation-dominated era in the early Universe. 

^Zf} , Causality dictates that the field lines form a tangled random network, and the string-like equations 

■ of motion, plus the assumption of perfect reconnection, inevitably lead to a self-similar solution for 
I ' the magnetic field power spectrum. We present the predicted form of the power spectrum, and dis- 
cuss insights gained from the string approximation, in particular the implications for the existence 

^ • or not of an inverse cascade. 
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^ ! I. INTRODUCTION 

> ; 

■ It has been observed that many galaxies and clusters of galaxies are endowed with a magnetic field with typical 
^y-j \ strength of order 10^^ G The origin of these large scale magnetic fields is unknown. In order for magnetic fields 

. to have this order of magnitude it is widely believed that an enormous amplification of an initial seed field must 
' have taken place. This amplification is usually explained by dynamo theory which can enhance the magnetic field 
^\ , exponentially However, dynamos cannot create a magnetic field, and so in order for them to act they require 
0^ ' a seed. At present it is not clear whether this seed field has its origin from some astrophysical mechanism after 
(-H ] recombination, during the epoch of galaxy formation and afterwards or whether the seed field is of primordial origin, 
Q_i' produced in the very early universe. In the latter scenario it is believed that the primordial magnetic field would have 
been frozen into the highly conductive plasma as the universe expanded and cooled. Because of the high conductivity 
diffusion would be small and magnetic fiux conserved. If a magnetic field was produced in the early universe and was 
present at the time of recombination it may have had a significant effect on many astrophysical processes including 
^ , the formation of galaxies and stars. 
ILV ■ There are several ways of obtaining limits on cosmological magnetic fields. Limits have been obtained by Faraday 
. I rotation measurements of intergalactic fields [HJJ] . Other constraints have been obtained through the consideration of 
' the effects of magnetic fields on primordial nucleosynthesis and on the distortion in the microwave background 

5h ] due to the presence of a cosmological magnetic field |^,^ . 
. . . 1 Even if a primordial magnetic field was to weak to be of astrophysical significance, it is still of principal interest 
to study cosmic magnetic fields today because they can provide direct and important information about the kind of 
physics that must have taken place in the early universe. There have been quite a few mechanisms proposed for ways 
of producing magnetic fields in the early universe. We will not discuss them here but refer to [p^ for a brief review. 

In this work we will not be concerned with any particular model for the generation of primordial magnetic fields. 
Instead we will focus on the universial problem of how a primordial field, whatever its origin, will develop as the 
universe evolved. In order to do this one needs to consider magnetohydrodynamics (MHD) in an expanding universe. 
Doing full numerical rclativistic MHD simulations of the physics of the early universe is hard and requires extensive 
computer memory and time. 

Greater dynamical range can be obtained by resorting to approximate methods. The cascade model of Ref. |l^] 
is one such method, which is thought to reproduce well the flow of energy between wavenumbers of the full MHD 
equations, at the cost of a severe truncation in the number of degrees of freedom. It was found in JTl] ! that energy was 
transferred from small to large scales in an inverse cascade, and that the correlation length of the initially random field 
increased with what looked like a small power of (conformal) time. This is pleasing if one wants to derive the galactic 
dynamo seed field from a primordial process, for general arguments of causality and energy conservation indicate that 
such an inverse cascade is actually necessary pO| ]. 
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In this work we will be using a string model approach to relativistic MHD to study the evolution of cosmic magnetic 
fields. The connection between MHD and string dynamics have previously been studied by Semenov |]l2| and Olesen 
psf . Our approach is simular to that of Semenov [|2j but more general since we do not assume that there is a 
conserved particle number density. We take essentially the opposite direction of Olesen |^ , in that we derive string 
equations from MHD and not the other way around. 

Once we have reduced the MHD equations to a string model, the results can be understood in terms of the 
coarsening dynamics of cosmic string networks [ p^HTs] ] . The rate of increase of the network scale length ^ is given by 
the characteristic velocity of waves on the string, in this case the Alfven velocity, which decreases as the magnetic field 
decreases in strength. The string approach indicates that this decrease in strength is primarily due to reconnection 
on small scales: small flux loops are continually created, transferring energy away from the network of infinitely long 
flux lines. The transfer of energy from the large-scale field happens in a self-similar manner: the magnetic field power 
spectrum can be displayed as 

|Bkp (X t-"P(A:C(t)), (1) 

where k is wavenumber and r conformal time. A powerful scaling argument due to Olesen [l^ shows, in the limit 
of ideal MHD, that <^ c>c r^/'^"^^^ where the initial power spectrum behaves as /c" at low k. Causality dictates that 

> 2 ^ (and not n > as one of us and another author has stated). As we violate the ideal condition by 
allowing reconnection, it is not clear that this is the correct power of t. 

This scaling law is our main result. We see no sign of a true inverse cascade, in the sense that power is not 
transferred from small to large scales. If anything, the transfer is from large to small, and it is only because energy 
is being lost faster from small scales that we see an increase in the scale length ^. 



II. RELATIVISTIC MHD AND STRINGS 



In this work we will concentrate on the ideal limit of MHD. This means that we neglect any viscous effects and treat 
matter as a perfect fluid. This is a good approximation at sufficiently large scales. Duringthe radiation dominated 
era, which we are mainly concerned with here, the universe was a very good conductor |l8||l9| ]. We therefore consider 
the CT ^ oo limit of MHD where magnetic diffusion can be ignored and the magnetic field can be considered to be 
frozen into the plasma and thus conserving magnetic flux. 

The starting point for ideal relativistic MHD is the energy-momentun tensor 

r/'^r/" 1 1 

47r 4 

consisting of the ideal fluid part and electromagnetic part of the energy-momentum tensor. Here p is the fluid 
pressure, pc^ is the energy density of the fluid, [/'^ is the four-velocity of the fluid satisfying the normalisation 
condition U^U^ = (? and F^^ is the electromagnetic field tensor. 
The evolution equations for the system are given by 

T^^^ = (3) 
47r 

i^^''.^ = — J'' (4) 
c 

*F^'',, ^ (5) 

aF^^U^ - rU^U^" (6) 

where equation (^) expresses covariant energy-momentum conservation, equations (jj) and (||) are Maxwell's equations 
with being the four-current density. In equation (||) *F^^ is the dual field tensor defined through the relation 



FflV — n^flV^pF^ (7) 



where e^y-yp is the Levi-Cevita symbol. Equation (|^) is the relativistic version of Ohm's law where a is the conductivity 
of the fluid, measured in the fluid rest frame. 

We now repeat the derivation in Subramanian and Barrow pc| ] to show that the evolution equations are conformally 
invariant and the evolution can therefore be transformed from the expanding universe to a flat (Minkowski) spacetime. 
In so doing we obtain an equivalent set of equations which are easier to handle. Two metrics g^i, and g^^, are said to 
be conformally related to each other if g^^ — fl^g^i, where 51 is a non-zero differentiable function. 
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The flat Robertson- Walker line element has the form 



ds^ = df -a^{t)d-x^ (8) 

where t is the comoving proper time and a{t) is the scale factor. 

This metric describes a isotropic and homogeneous universe with zero curvature. The appearence of a hypothetical 
primordial magnetic field of some strength need not violate the assumption of isotropy and homogeniety because 
although the presence of a magnetic field will locally generate bulk motions in the fluid, if we look at sufficiently 
large scales isotropy and homogeniety will be regained. At large scales the magnetic field, whatever its origin, can be 
considered as essentially random since the correlation length of the field is bounded from above by not exceeding the 
causal horizon. This justifies the use of the Robertson- Walker metric. 

We introduce conformal time r defined by dr = a^^dt so that equation becomes 

ds^ ^a^{T){dT^ ~ dyi^) (9) 

and hence rj^i, = g^i, = ^'^g^^i, with = a^^(T). 

We note that under conformal transformations the ideal energy-momentum tensor T'"' transformation as T^'' = 
a~^T^'^ . That this is so can be seen directly from the definition of the energy- momentum tensor 

T'"' = iV^Lrnatter) (10) 

The new scaled fields (denoted by tilde) obey ordinary energy-momentum conservation. To see this we note that 
the ideal energy-momentum tensor is traceless, T = T'^^ = 0, provided the perfect fluid has the equation of state 
p = \pc^ ■ For most of the period before decoupling, the early universe was radiation dominated and one can use the 
above equation of state. We have 



Using 



we get 



But since T = then 



and hence 



T>:,^^-^{5^Jl + 5^X-g„X) and TJ^^ = 4^(5^) (12) 



^ = % + 2-T'^° - -5^T + 4-r''° (13) 
a a a 



(a^^f ''0) ^ + 6-(a-6T'^°) = {f^"" Jfl-^ = (14) 
a 



f""^^ ^ (15) 

This means that under conformal transformations our original variables will transform to a set of new scaled variable 
satisfying the following relations: p — a^^p, p = a^^p, = a^^U^, = a^'^J^, F^^ = a^^F^^ . 
Now consider Maxwell's equation (^). Since F^" is anti symmetric the left hand side simplifies to 

F^"',p F^-'^p + T^^F'"' = a-^F^"'^p (16) 



So the equation for the scaled fields becomes 



F"" , = — > (17) 
c 



For the four-velocity we have 



= g^pUP = a'g^pia-WP) = all, (18) 
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Hence 



aF^'^U^ = aa-^Fi"'U^ = a'^iJ^ + J'^U^U^) (19) 

So Ohm's law remains invariant under conformal transformations if we define the scaled conductivity through a = aa~^ 
So we arrive at the fundamental equations of relativistic MHD, 

f""^,^ = (20) 

F^"'^ = —J'' (21) 
' c 

*F^"'^^ = (22) 

aF>'"U^ = jf" - rU^U^'/c^ (23) 

Prom here on we drop the tilde, on the understanding that we mean scaled fields. 

Wc will now introduce a new set of coordinates which will enable us to write the MHD equations as non-linear 
string equations. We define a magnetic four-vector through the relation 

6" = *F^.U''/c (24) 

We also define new coordinates x' = {rj, a, t/j, () such that 77 are coordinate lines of fluid elements and cr arc coordinate 
lines of magnetic flux. Hence = {cq, 0, 0, 0) and b'f" = (0, P, 0, 0) satisfying U'"^ = = and b'^ = 6^ = -B'^ 
respectively. Thus we have the metric tensor in the new coordinates 

g'^, = diag{l/q^,-B''/p\-hAB) 
where A,B = 2, 3. The new coordinate vectors are 

^ = — (25) 
or] q 



^ = 1 

Since the introduction of these coordinates relies on the frozen in property of the plasma we will refer to them as 
frozen-in coordinates. In the frozen-in coordinate system we can trace the trajectories of fluid elements by simply 
varying the value of our time coordinate 77 and keeping the values of the other coordinated fixed. Similarly, wc can 
trace the magnetic field lines in the frozen-in system by varying the value of a and keeping the other three coordinates 
at fixed values. 

The analysis of the MHD equations is usually performed in terms of the magnetic field and velocity distributions. 
However, in the description of MHD phenomena, the concept of a magnetic flux tube is often introduced. The reason 
is that it is sometimes convenient, and we gain a better physical insight, if we base the description on this concept 
rather on the magnetic field and velocity distributions. A magnetic flux tube is defined as the volume V enclosed by 
a closed surface E which is everywhere parallel to the ambient magnetic field vector, and two cross-sectional surface 
areas Sa and Sb at either end. The flux tube therefore consists of a bundle of magnetic field lines which enter and 
exit the volume through the end of surfaces Sa and S b ■ 

The Cartesian coordinate system, traditionally employed in MHD does not really lend itself to an analysis of the 
magnetic flux tube bahaviour. The frozen-in coordinate system on the other hand, docs. The frozen-in coordinates 
provide a coordinate system co-moving with a flux tube, and is therefore the more natural choice for the mathematical 
analysis of flux tube behaviour. 

We now consider the equation for magnetic evolution 

*Fi^''.^ = (27) 

Since *F^" is antisymmetric the divergence is given by 

1 



^^{^/^''F'"') = (28) 
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Using the fact that we can express the dual field tensor as 
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(29) 



we have 



which gives 



1 / 6^ 

But first square bracket is the Lie derivative and so vanishes, giving 



Thus we have 



-rdAV^PUn dAV^qb'') = 



(30) 



(31) 



(32) 



and 



Hence in our comoving frame we get 



|^(V^/39) = ^(Bdet(/i^B))=0 



(33) 



(34) 



(35) 



d 



da 



And so we see that 



d 



^iV^M - ^{BdetihAB)) = 



da 



(36) 



Bdet{hAB) = F{^,C) 

where F{ip, () is an arbitrary function of ip, (. We therefore have the freedom to choose det{hAB) such that F{iIj, (^) = B 
where 5 is a constant and so 

hAB — -f^^AB 

D 

which means that we can write 



_ B 

We now study the equations of motion, starting from the energy-momentum tensor given in equation (^. Using 
the above expression, equation (|2^), for the dual field tensor and the connection between the field tensor and its dual, 
equation (^, we can write the energy- momentun tensor in the following form 



Ait 



47r 



r>2 



(37) 



From energy-momentum conservation we have 
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m ■ d) 



(38) 



Here P is the total pressure from both fluid and electromagnetic field. Note that the first and the third terms in this 
equation vanish. Writing this equation in our comoving reference frame we find 



!3q 



dr] \ V 



pc +P + 



q 1 d f P dxt' 

in ) c[3 drj ) At: da \q da 



'P.. = 



(39) 



Equation ( |3S| ) is the equation of motion in the frozen-in coordinates. In the frozen-in coordinate system the MHD 
equation of motion reduces to a set of non- linear string equations. The behaviour of a magnetic flux tube is therefore 
formally analogous to that of a non-linear string. The last term of the left hand side of equation ( |39| ) take account of 
inhomogeneity (i.e. pressure gradients) and it describes the coupling between neighbouring flux tubes whilst moving 
through a non-uniform plasma medium. 

To summarize, we have shown that the behaviour of a magnetic flux tube is formally analogous to that of a string and 
one can therefore model a magnetised plasma as a fluid composed of strings. We have relied heavily on the arguments 
of Semenov and Semenov and Berkinov Jl^ , with one improvement: we have dropped their assumption that there is 
a conserved particle number density, which is neither necessary nor generally applicable in the early Universe. Our 
derivation is also complementary to that of Olesen who starts from the relativistic string equations and shows 
that they can be interpreted as describing the motion of narrow flux tubes, providing the total pressure remains 
constant across the tube. 



III. APPROXIMATE STRING EQUATIONS 

Exploiting the freedom to change coordinates in the cr, r] subspace, we choose (3 such that (5 = qB and using this in 
the equation of motion (B3) we have 



1 d fcBdx^'\ 1 d f^dx''\ Aw 
Bdr]\j^~&q J ^ Bd^ \ ~d^ J ^ 



l^a^'^P,'^ - 0, (40) 



"A 

where we have deflned the relativistic Alfvcn velocity as 



cB/V4^ 

= 77?- (41) 



Rearranging equation ( pO| ) we can write it as 

1 a / 1 ^x^'\ d^x^' 1 d f cB\ dx>' 1 dB dxf" in 



g'^'P,. (42) 



va drj \VA drj J da"^ B drj \ va J drj B da da B'^q^ 

We now argue that for the particular situation we are interested in, it is justified to neglect the three terms on the 
right hand side of equation ( 4^) . The third term on the RHS of equation ( p2| ) can in general not be neglected because 
in many astrophysical situations pressure gradients are important. However, in the early universe the pressure from 
the magnetic field should be much smaller than the fluid radiation pressure and since the early universe had a very 
high degree of homogeniety it follows that the gradients of the total pressure are small. 

The second term on the RHS of equation (^) does not have have a definite sign and so if averaged over time it will 
be zero. Dropping this term is equivalent to replacing B by its root mean square value. 

Again using the fact that the fiuid radiation pressure in the early universe was much larger than the pressure from 
the magnetic field and remembering expression (^l|) for the relativistic Alfven velocity it is seen that the first term 
on the right side of equation (|4^) is indeed small since the ratio B /va is approximatly constant, thus justifying our 
decision to neglect it. 

We now rescale the time parameter 77 to f] through the relation 
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d c d 
dfj va{v) dif 

where we have called attention to the fact that the Alfven velocity may be time-dependent. Using the above mentioned 
approximations and our new time parameter 77 we are left with the following equation of motion 

^ 0. (44) 



We also recall that the four-dimensional orthogonality of the coordinates, and the coordinate choice enforced by 
/3 — qB, supplements this equation with constraints 

(45) 

where dot (prime) denotes derivative with respect to the timelike parameter fj (a). 

We note that the equation (44) and constraint (|5|) is of exactly the same form as the equation for a Nambu-Goto 
string in Minkowski spacetime expressed in the conformal gauge [Q^. With ideal Nambu-Goto strings in Minkowski 
spacetime one can further take x'^ = cfj, to obtain the system 

^x - x" = 0, (46) 
X • x' = (47) 
^jc^ + x'2 = 1 (48) 

This is not possible in general for MHD strings, for one can easily verify that 

V • B 

x-x' = c2j— -, (49) 

lx^-Hx^ = l!f^+(^)' 
(? q^ \'^\ c^B^ 



where 7 = (1 — v/c) ^1"^ is the usual relativistic gamma factor. However, as long as v • B = in the initial conditions, 
the first constraint is preserved by the evolution. Furthermore, we can use our remaining coordinate freedom to define 
9 by 

9' =7'^, (51) 



in which case we really do reproduce the Nambu-Goto equations. We should bear in mind however that we have made 
several approximations on the way, which are worth reiterating. 

1. We have neglected pressure gradients. 

2. We have made a kind of mean-field approximation in replacing the magnetic field by its root mean square value. 

3. We have neglected the effect of the time-dependence of B/va- 

We have argued that all this approximations are reasonable in the context of the early Universe, and find that the 
Nambu-Goto equations can be used to approximate a class of MHD velocity and magnetic field configurations which 
satisfy v • B = 0. They may also be a reasonable approximation to other configurations, provided we understand that 
the constraints are satisfied in an average sense. 



IV. ALGORITHM, SIMULATIONS AND RESULTS 

The equation of motion (^6|) can be evolved using the Smith- Vilenkin algorithm . The Smith- Vilenkin algorithm 
provides an exact discrete evolution for a string network defined on a face-centered cubic lattice and the evolution 
equations are 
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x{f] + S,(t) = ^ [x(r;, (T + S)+ x{r], a - S)] + v(r;, a)6 (52) 

v{f] + S,cr) = ^ [v{i], (j + d) + v{7], a -6))] + ^ [x{f], a + 2S) - 2x(r;, a) + x{r], a - 25)] (53) 

By using a leap-frog method, updating the positions and velocities at alternate timesteps, the above equations allow 
us to calculate the exact future evolution of x and v from some appropriately choosen initial conditions. 

Initial string configurations are generated by a method due to Vachaspati and Vilcnkin . They considered string 
formation in a global U{1)- model. The J7(l)-manifold is discretizcd by allowing the phase to take only three possible 
values. These values are then placed randomly on the sites of the cubic lattice. As we go around the face of a cube in 
real space, the phase descibes a certain trajectory in the group space. A string passes through the face of a cube if that 
trajectory has a non-zero winding number. With this method the string segments join up to form either closed loops 
or else open strings which intersect the boundaries of the cubic lattice. The shape of the strings will be Brownian 
with step size corresponding to the cell size of the lattice. 

We have seen how relativistic MHD under the approximations discussed above can be described in terms of magnetic 
flux tubes satisfying the Nambu-Goto equations of motion. In order to represent the continuous distribution of 
magnetic flux by a network of string, we gather together a flux <& into ideal Nambu-Goto strings at positions X*(^, a), 
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B\t, x) = * I da ^^'l^yh ^ (x - X(fy, a)) (54) 

A few words need to be said about reconnection. In real fluids, magnetic flux tubes interact and reconnections 
take place when the field lines cross each other. In ideal MHD there are no dissipative or viscous effects. Physical 
reconnection between field lines cannot take place without resistive effects and therefore the topology of the magnetic 
field lines is frozen in the fluid and does not change with time. 

Reconnection, which is a local process, is difficult to describe and the details are not well understood. The efficiency 
of reconnection is not known and neither is how this process depend on the local geometry involved, like the relative 
inclination of the flux tubes. To answer this question, one has to go to numerical solutions of the underlying theory, 
that is MHD. 

In order to allow for reconnection to take place in our simulations we will take a more simple approach to re- 
connection. In our model strings are allowed to reconnect only if they pass through the same lattice point. When 
two strings meet, they intercommute with a probability P and in this work we put P = 1. Reconnection here take 
place instantaneously between the discrete evolution time steps. This is physically reasonable since the reconnection 
timescale is small compared with the evolution of the system. 

We use a simple estimate for the characteristic length scale ^ which is defined by 

e = V/L, (55) 

where V is the volume of the lattice and L is the total length of string in the box, excluding small loops. In the 
Smith- Vilenkin algorithm it is in fact possible to have a loop of zero spatial extent, occupying just one lattice point, 
which travels at the wave velocity. If such a loop is formed it does not contribute to the magnetic field, and we do 
not count it in the calculation of the total length. It is a well-known feature of Nambu-Goto simulations that nearly 
all the string ends up in this kind of loop In the MHD context we should probably not regard these loops 

as representing magnetic field lines at all: the energy is probably being dissipated in a very small-scale reconnection 
process. In any case, the fact that loops are generically produced so small underlines the fact that strings are very 
efficient at transferring power from large to small scales. 

We have studied simulations of the evolution of the magnetic field on lattices with sizes (1286)^ and (256(5)^, using 
a code originally developed by Sakellariadou which implements both the Vachaspati- Vilenkin algorithm for the 
initial conditions and the Smith- Vilenkin alorithm for the evolution. Periodic boundary conditions were used in all 
simulations and the evolution time were resticted to less than half the box size, since after this time causal influences 
have had time to propagate around the box. 

Having a representation for the magnetic field in real space we use a three dimensional Fast Fourier Transform 
algorithm to get a Fourier mode representation. The power spectrum for the magnetic field |Bkp can then be 
calculated at every time step by averaging over the amplitudes of all Fourier modes with a wave number between k 
and k + 2Tr/S. 

Previously it has been shown that networks of cosmic strings, modelled as Nambu- Goto strings in Minkowski 
spacetime tend towords a scaling regime M . This means that the characteristic length scale of the network grows as 
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^ ct where t is Minkowski time. The characteristic length scale for a primordial magnetic field ^ does not grow with 
the horizon. Instead we expect the magnetic field to grow as ~ C77. 

It is interesting to see if the magnetic field power spectra show scaling behaviour. In order to investigate this we 
express the power spectrum in terms of a scaling function P which is defined through the following expression 

\B{kM'=V^^^. (56) 

Here V is the volume of the box and its appearance is just a normalisation convention. This form for the scaling 
function ensures that that the fluctuations obey the scaling law 

which is consistent with our picture of a coherent flux $ in a region of size ^. 

Figure 1 shows the chactcristic length scale of the magnetic field versus cfj for a typical ensemble. The measured 
scaling function P(fc<^) is displayed in Figure 2. The data was taken from 12 runs on a (2566)^ lattice, for values of 
cfj between 40 — 60. It is seen that the power spectrum of the magnetic field does reach a scaling regime. This means 
that the evolution of the network will be self-simular with respect to cfj. 

It is clear that the network does exhibit the property of scaling, with — oix, ■ cfj. The scaling amplitude can be 
obtained by looking at the ratio £^/fj towards the end of the simulation and it is roughly a;. ~ 0.3 (see for a more 
accurate determination). Thus we can write 

S, — ■ cfj — Xt, J VA{r])dri. (58) 

Unfortunately, we do not yet know how va depends on the time parameter rj or the conformal time r. The Alfven 
velocity va depends on _B, which in turn depends on $ and ^ through (^6|). All we can infer from the information at 
hand is a consistency relation: if ^ oc 77'' and $ oc yy'', then 3r = s + 1. 

The extra information we need comes from the covariance of ideal MHD under the scale transformation |jl^ 

x->Zx, T~^l^-^T, v^l^v, B^l'^B, (59) 

where h is arbitrary. One can show that under this transformation, 

y-i|Bk|^ ^l^+'^V-i|Bk|^ (60) 

If we define a function A(fc, t) by 

l^-i|Bkr = T(3+^)/(^"'^)A(k,r), (61) 
then we see that under the same transformation 

A(fc, r) ^ A{k/l, tI^-'') = A{k, t), (62) 

from which we immediately infer that 

A(fc,T) = A(/ct1/(i-'*)). (63) 

Furthermore, if |Bkp behaves as fc" as A: ^ 0, we have 

|Bkr cxr(^+'^+")/(i-'^)k" (64) 

in that limit. It is often assumed that the large-scale power is not affected by small-scale processes |16 1^, in which 
case h = —n — 3, and we find the scaling laws (derived by the same authors) 

^cxt2/("+5), $oct(i-")/("+5). (65) 

In the early universe individual particles move with relativistic velocities. However, we expect that the bulk velocity 
of the fluid v to be non-relativistic. Hence 

[/" = qc^ ^ c, (66) 



and given ( |51| ) we see that, on average, 

rj^^{k^T/vA, (67) 
where (x^) is the mean square string velocity. Simulations give this to be 0.36 123. 
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V. CONCLUSIONS 



We have seen how the relativistic MHD equations, with a few reasonable assumptions, may be recast as string-hke 
equations for the motion of the flux hnes. This allows us enormous gains in dynamic range in the simulation of a 
random magnetic field in the early Universe, without being forced to the ideal limit, for we incorporate diffusivity by 
allowing reconnections between the strings. 

The result is that we can understand the evolution of magnetic fields in terms of the evolution of a network of 
strings, and we find that the power spectrum quickly evolves to a self-similar or scale-invariant form, with scale length 
f increasing in time. What this power law is we are unable to say: ideal MHD predicts f (x t^/*^"+^\ where n is the 
low k exponent of the power spectrum, but as we have departed from ideality by allowing reconnection, we cannot 
make a prediction. 

The increase in scale length comes about by the strings straightening at the Alfven velocity, while forming very 
small loops which can dissipate energy quickly. This is the new feature that the string formulation brings to light: 
strings transfer energy from large to small scales in an extremely efficient manner. Thus, although the scale length 
increases, it is because power is preferentially lost from small scales. Whether it is fair to call this an inverse cascade 
is a matter of terminology. What is clear is that the dynamics predicted by the string model of MHD is certainly not 
of the right kind to produce seeds for the galactic dynamo from magnetic fields created in the very early Universe. 

There are of course many places where this line of argument is vulnerable. The model makes approximations which 
we have tried to highlight. Furthermore, our string simulations use special string configurations to make gains in 
computational efficiency: the strings lie on a cubic lattice to start with, and one may be suspicious that this may 
introduce some artifice into the dynamics. However, the propensity of a string network to scale is firmly believed, 
so we are confident that the magnetic field power spectrum will also scale. What is probably not well approximated 
is the actual form of the power spectrum, which betrays the particularly string-like feature of a tail, due to the 
fact that all the fiux is held to be concentrated in a narrow tube. Furthermore, the string model may be deficient 
in its description of helicity, which is known to be extremely important in the development of true inverse cascades 
pSj-pTf. The helicity is represented by the linking number of the strings, but we are not able to incorporate a local 
contribution induced, for example, by twisted tubes of flux. It may well be that we are missing some very important 
dynamics here. We clearly need to check our results against a non-ideal MHD code, to see if the predicted self-similar 
dynamics emerges, and also to find the correct power law for the scale length. This project is currently in hand. 
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